Smoothing of all-sky survey map with Fisher-von Mises function 
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Convolution of all-sky survey data with a smoothing function is crucial in calculating the smooth 
surface brightness of the sky survey data. The convolution is usually performed using the spheri- 
cal version of the convolution theorem. However, a Gaussian function, applicable only in flat-sky 
approximation, has been usually adopted as a smoothing kernel. In this paper, we present an 
exact analytic solution of the spherical harmonic transformation of a Fisher-von Mises function, 
the mathematically version of a Gaussian function in spherical space. We also obtain the approx- 
imate solutions exp[— 1(1 + l)/2«]. The exact and approximate solutions may be useful when an 
astrophysical survey map is convolved with a smoothing function of k > 1. 
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I. INTRODUCTION 

In order to compare the results of theoretical calcula- 
tions with an astrophysical observation, it is necessary 
to convolve the theoretical calculations with a measure 
of the angular response of the detector and to calculate 
the quantities which correspond to those being observed. 
The angular response is often called the "beam function" . 

In particular, knowledge of the beam profile and its 
spherical harmonic transformation is of critical impor- 
tance for interpreting data from cosmic microwave back- 
ground (CMB) experiments. A common simplifying as- 
sumption in the spherical data analysis is to take the 
experimental beam response to have a circular Gaussian 
profile as in the flat-sky approximation [l], Q ■ More re- 
cently, methods to perform the convolution of a CMB 
map with a general asymmetric beam, an elliptical Gaus- 
sian beam profile have been also investigated 0, 0] ■ 

It is often required to convolve spherical data, such as 
all-sky survey maps, with a smoothing function. A cir- 
cular Gaussian beam function has been also used to take 
the angular position errors into account in the procedure 
of computing the angular power spectrum of BASTSE 
3B gamma-ray bursts § • Smoothing of two-dimensional 
spherical data is also a procedure routinely used in many 
fields of sky survey data analysis [f|. 

A Gaussian function in flat-sky approximation is ap- 
plicable only when the beam width a <C 1 . The spherical 
version of a Gaussian function is known as a Fisher-von 
Mises function by mathematicians @. Fig. 1 shows sub- 
stantial differences between a Gaussian function and a 
Fisher-von Mises function when the beam width a < 1. 
Thus, it is worth to investigate the convolution of spheri- 
cal data with a Fisher-von Mises distribution function. In 
this paper, we investigate the spherical harmonic trans- 
formation of a Fisher-von Mises function, which is re- 
quired for the convolution of a spherical data with a beam 



width a < 1. We found the exact solution, recursive re- 
lation, and generating function of the spherical harmonic 
coefficients. 

The FIMS (FUV IMaging Spectrograph) instrument 
was launched onboard the first Korean Science and 
Technology SATellite STSAT-1 on 27 September 2003 
[USES]- Its spatial resolution is about 10' — 30' due to 
the optical properties of the FIMS instrument and errors 
in attitude knowledge provide by spacecraft. The first 
far-ultraviolet (FUV) sky maps are in preparation to be 
published. Analysis of the FUV sky survey data obtained 
with FIMS motivated the present work. 



II. SPHERICAL HARMONIC CONVOLUTION 

We describe the spherical harmonic convolution the- 
orem in order to define naming conventions. Let the 
continuous function A{6, <ft) represent a map of the sky, 
and let B(9) be the (azimuthally symmetric) smoothing 
kernel. If we expand a sky map A(9, </>) in spherical har- 
monics as 



A(0, 



E E Al m Yl 



1=0 m=-l 



the coefficients Ai m arc given by 



A, 



Y l * m (9, ( / ) )A(0^)dn. 



(1) 



(2) 



Similarly, let us expand the beam function in Legendre 
polynomials as 



2=0 



B{9) = Ya/^P-Bu.Y,;^.- 
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Here, we should note the different convention for the co- 
efficients Bi. Because of azimuthal symmetry of B, the 
Bi m coefficients vanish for m ^ so there is no sum of m, 
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FIG. 1: Comparison of Gaussian and Fisher- von Mises func- 
tions with beam widths a — 0.1,0.3,0.5, and 1.0. Solid and 
dotted lines denote Fisher-von Mises and Gaussian functions, 
respectively. 



and we abbreviate Bio with B\. Inserting the definition 
of Y w gives 



B(9) = £('2*±T) 



4-7T 



The coefficients Bi are, then, given by 



Bi= I B(6)Pi(cos9) sin 6d6. 
la 



(3) 



(4) 



The convolution of A with B is defined by 

H{9, cj>) = A{9, <j)) * B{6) = [ dn'A{6', 4>')B(i), 



where 7 is the angle between (9, <fi) and (6 1 , <f>'). If we de- 
note the spherical harmonic transform of H(9, <fi) by Hi m 
following the definition Eq. O, we have the spherical 
harmonic convolution theorem [5|, 



Hi, 



BlAfr, 



(5) 



III. SPHERICAL HARMONIC 
TRANSFORMATION 

A Fisher-von Mises function 0] is defined by 
Kexp(ftcos#) 



B{6) 



47T sinh(K) 



(6) 



characterized by the concentration parameter k. The 
larger the value of k the more the function is concentrated 



towards the direction 9 = 0. The concentration parame- 
ter is related to the usual beam width by k = 1/a 2 . This 
is often considered by mathematicians to be the spherical 
version of a Gaussian distribution, and it reduces to 



B(9) = 



exp(-6» 2 /2cr 2 ) 



(7) 



2vrcr 2 

when (r < 1 or k > 1. This is the smooth kernel which 
has been commonly used for the convolution of all-sky 
survey maps in literatures [3, 0, ■ The Fisher-von Mises 
function has the advantage that it is correctly normalized 
(its integral over the sphere is unity) for arbitrarily large 
beam width a, which is not the case for the plane Gaus- 
sian of Eq. 0. 

Substituting x = cos 9 into Eq. §5§ , we obtain 



B, = 



+1 d xPl ( X ) KeXp{Kx) 
-1 



2 sinh k 

Integrating this by parts, we obtain 



(8) 



B, = 



9 \e K Pl{l)-e- K Pi{-l)- [ + dxe™P{ 

2 Sinn /c |_ J_ 1 



(x) 



(9) 



Using the recursion relation P[ +1 (x) — P/_ 1 (a;) = (21 + 
t)Pl(x) and the properties P/(l) = 1 and P/(-l) = (-1)' 
[111 ] , we have the recursion relation for the spherical har- 
monic coefficients P/, 



B, 



1+1 



-(21 + 1)Bi/k. 



(10) 



The first two coefficients, Bo and B\, can be easily ob- 
tained from Eq. §§§ as 

Bo = 1, and 

B\ = cosh« — 1/k. 

Using theses, we have also the following recursion rela- 
tion. 



(11) 



Now, we obtain an explicit formula for the coefficients B\ 
by continuing the integration by parts. 



B, = 



1 



1 

y 

2 sinh k *-~> 

,(n) 



3 KX P/ n) (l) 



p(n) ( _ 1} 



Here, P z = d n P/dx n . Using the Rodrigue's formula 

El 



1 d l 



we have P/ n } (-l) = (-l)'+ rl P/ n) (l). Thus the formula 
for the coefficients Bi is given by 

*-ia E ;X:(~)"[T--(- i n'f , (i) 

n=0 v 7 
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In order to calculate P[ n \l), we substitute x 2 — 1 
2(x — 1)[1 + (x — l)/2] into the equation and obtain 



1 cz'+ n 

Z!tfe I +™ 



(x-iy i + 



X - 1 



1 1 1 

_ v - - ( 

cfe'+ n ^ 2 m (Z -to)!to' v 



x — 1 



>A 1 1 (Z + m)! 
^ 2 m (I- m)\m\ (m - n)\ 



(x-iy 



Thus, we have 



p(")(!) - 1 ( l + n ) [ 



(12) 



2" (Z - n)!n!" 

Further, we can rewrite the formula when Z > is odd 
Bl = coth. E f 1 ) ^C 1 ) 

m=0 ^ ' 
(2-l)/2 . . 2m+l 



- E (i) if—m 

rn—O 
n=0 v 

+(coth«-l) E (£) ^ (2 " l) (l) 



P w (1) 



(Z-l)/2 



m— 



When Z > is even 

2/2 / , x 2n, 



Bl = e (i) ^ (2m) (D 

m=0 ^ ' 



V 2 -! / -, \ 2m+l 

-coth.E (~) ^ m+1) (l) 



71=0 



1/2-1 



-i - - / 1 \ 2m+l 

(coth«-l)E ^/^(l 



When k y§> 1 , coth k w 1 and we have 



B < - Ef-;)V n) (D 



n=0 



n=0 



(i + n)! 
k/ (/ - n)\n\2 r < 



Since (/ + n)\/{l - n)\ « [Z(Z + 1)]" for Z > n, we have 



5/ 



1\" W + 1)] T 



i!2* 



exp 



Z(Z + 1)" 
2k 



(13) 
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FIG. 2: Comparison of the spherical harmonic coefficients for 
Gaussian and Fisher-von Mises functions with beam widths 
(a) a = 0.1, (b) 0.3, (c) 0.5, and (d) 1.0. Solid and dashed 
lines denote the coefficients calculated for Fisher-von Mises 
and Gaussian functions, respectively. Dotted line represents 
the approximation exp[— 1(1 + 1)/2k]. 



This is the same approximation as calculated for a plane 
Gaussian function [H, H, [H . 

In Fig. 2, the spherical harmonics coefficients for 
Fisher-von Mises functions are compared to those for 
plane Gaussian functions. The coefficients for Gaus- 
sian functions were calculated numerically. The approx- 
imation Eq. (|13|) are also shown for comparison. It is 
clear that the approximation represents the coefficients 
for Fisher-von Mises functions very well when a < 1. On 
the other hand, the Eq. (fTBf approximates the coefficients 
of Gaussian function only when uCl. 



IV. GENERATING FUNCTION 



The recursion relation Eq. (|10[) is a difference equation, 
thus we can obtain a differential equation equivalent to 
the difference equation [12j , by defining that 



F(x) = E 



1=0 



Bix L 



(14) 



To obtain the differential equation for F{x) we multiply 
Eq. (PH by x l ~ l /{l - 1)! and sum from I = 1 to oo. The 
first term gives 



E* 

i=i 



.i-i 



i+i 



(Z-l) 



dx 2 

d 2 F{x) 
dx 2 



[F{x) -B - B lX ] 



The second term gives 



J-i 



E ^-1(731) 



F(x) 



i=i 
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The third term gives 



We may expand the generating function as, when k<1, 



Combining these three results gives the differential equa- 
tion for F(x) 



(1 + 2x/k)F" + 3F'/k -F = 0. 
This equation is simplified by substituting y 



(15) 
(1 + 



k 2 dy 2 
Thus, we obatin 



y d 2 F 2 dF 
n 2 dy 



yF = 0. 



ht* {vF) - yF=Q - 

The general solution of the above equation is given by 

yF{x) = de Ky + C 2 e~ Ky . 

Using the initial conditions, F(0) = B = 1 and -F"(0) = 
B\ = cothK — 1/k, we obtain the generating function for 
the spherical coefficients Bi, 



F(x) 



sinh (fty/l + 2x/kJ 



(sinh k) \fl + 2x/k 
The formula for Bi is, then, given by 

d l F{x) 



B, 



dx l 



(16) 



(17) 



F(x) = - 



i ^ n 2n+1 (i + 2x/ K y 



sinh« ^— ' (2n + 1)! 

n=0 K ' 



Thus, we obtain the coefficients B{ for k <C 1 



By = 



(2/«) 



EK 1 " 
(2n+ 1)!' 



sinh/c ^-J (2n+ 1)! (n - Z)! 
2'Z! 



(2Z + 1)! 



SUMMARY 



A recursion formula for the spherical harmonic trans- 
formation coefficients of a Fisher-von Mises function, 
which is the spherical version of a Gaussian function, was 
obtained. The formula may be used for the convolution 
of all-sky survey maps, such as CMB maps and FIMS 
FUV sky maps, to obtain smoothed sky maps. 

We also investigated the mathematical properties of 
the spherical harmonic coefficients and found that the 
coefficients approaches exp[— 1(1 + 1)/2k], and cx k 1 when 
k > 1 and k<1, respectivley. 
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